Field-theoretical renormalization group for a flat two-dimensional Fermi surface 
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We implement an explicit two-loop calculation of the coupling functions and the self-energy of 
interacting fermions with a two-dimensional flat Fermi surface in the framework of the field theo- 
retical renormalization group (RG) approach. Throughout the calculation both the Fermi surface 
and the Fermi velocity are assumed to be fixed and unaffected by interactions. We show that in two 
dimensions, in a weak coupling regime, there is no significant change in the RG flow compared to 
the well-known one-loop results available in the literature. However, if we extrapolate the flow to a 
moderate coupling regime there are interesting new features associated with an anisotropic suppres- 
sion of the quasiparticle weight Z along the Fermi surface, and the vanishing of the renormalized 
coupling functions for several choices of the external momenta. 
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I. INTRODUCTION 

A better understanding of the physical properties of 
highly interacting electrons in two spatial dimensions 
(2D) is believed to be central for high-Tc superconductiv- 
ity. Soon after the discovery of the high-Tc superconduc- 
tors, Andersoni suggested that a strongly interacting 2D 
electron gas should resemble a ID Luttinger liquid state. 
This question remains unresolved to this date. Thanks 
to the high precision angular resolved photoemission ex- 
periments performed in a variety of materials 2 , we know, 
at present, important facts concerning the Fermi surface 
(FS) of the cuprates. The FS of underdoped and opti- 
mally doped Bi2212 and YBaCuO compounds contains 
both flat and curved sectors 3 . As a result, they are nearly 
perfectly nested along certain directions in momentum 
space. As is well-known, whenever there is a flat FS, the 
corresponding one-electron dispersion law resembles a ID 
dispersion. 

The cuprates are Mott insulators at half-filling which 
become metallic at very low doping^. At half-filling, 
Hubbard-like models have a square shape FS imposed 
by electron-hole symmetry. The FS changes as we vary 
the filling factor and, as soon as it is lightly doped, it ac- 
quires non-zero curvature sectors in k-space. However, in 
the immediate vicinity of half-filling, there are, at most, 
isolated curved spots in momentum space. Consequently, 
in a first analysis, one may neglect their presence alto- 
gether. Following this scheme, several workers investi- 
gated the properties of a 2D electron gas in the presence 
of a totally flat FS 5 i 6 i 7 i 8 i 9 . In their approaches, the FS 
is always kept fixed and never deviates from its original 
flat form. Besides that, their results conflict with each 
other. Conventional perturbation theory calculations , 
parquet method results^, as well as one-loop perturba- 
tive RG calculations? indicate that, for repulsive inter- 
actions, there is never a Luttinger liquid state in 2D. 
For the repulsive Hubbard interaction, the one-loop cal- 
culations indicate that the antiferromagnetic spin den- 



sity wave is the dominant instability in two dimensions. 
In contrast, applying bosonization methods, Luther was 
able to map the square FS onto two sets of perpendicular 
chains 8 . As a result of that, the corresponding electron 
correlation functions become sums of power law terms 
with exponents only differing in form from those of a 
Luttinger liquid 9 . Very recently, Rivasseau and collab- 
orators using a mathematically rigorous renormalization 
group analysis^ proved that, at finite temperatures, the 
half-filled Hubbard model in two dimensions with a per- 
fectly square Fermi surface is indeed a Luttinger liquid 
apart from logarithmic corrections. That result adds new 
elements to the possible existence of a Luttinger liquid 
phase in two dimensions. 

We report in this work a two-loop field theoretical RG 
calculation for the electron gas in the presence of the 
same flat FS model as used by Zheleznyak et al. 6 . Ex- 
perimentally, such a FS was observed quite recently in 
an ARPES experiment performed in La,2- x Sr x CuOi-$ 
(LSCO) thin epitaxial film under strainii. On the theo- 
retical side, all one-loop RG calculations presented so far 
for a flat 2D Fermi surface indicate a flow to a strong cou- 
pling regime. One should therefore check if a two-loop 
calculation changes that result or not. Besides, there 
have been recent works addressing the renormalization 
of the quasiparticle weight Z due to interactions^ 2 *^. 
Those works are not totally consistent since they ignore 
the feedback produced by Z in the RG equation for the 
coupling functions. We show that, as long as we restrict 
ourselves to the weak coupling regime, Z is only altered 
slightly from its unity value. As a result, in that limit 
the RG coupling flows are not severely modified by Z . 
In contrast, if we consider moderate couplings, there are 
substantial changes produced in the RG flow as opposed 
to the one-loop case, and Z indeed goes to zero in this 
region of coupling space. 

In our analysis, we assume that the FS and the Fermi 
velocity vf are not renormalized by the interactions. In a 
later work, we show how the renormalization of both vf 



2 



and the FS may affect our results. The novel aspect of our 
work, apart from using a field theoretical RG approach 
to deal with a 2D fermionic problem, is the fact that we 
explicitly calculate all the two-loop diagrams for the self- 
energy and the renormalized coupling functions, taking 
full account of the variation of the momenta of the coun- 
terterms along the FS. As a result, we are able to calcu- 
late the nonleading logarithmic contributions originated 
by the so-called nonparquet diagrams and the self-energy 
corrections which lead to a momentum dependent quasi- 
particle weight Z(p). To our knowledge, it is the first 
time that such a systematic and detailed two-loop calcu- 
lation in two dimensions is presented in the context of 
the fermionic RG approach. For simplicity, and to make 
easier the comparison with other RG results, we use the 
so-called "g-ology" notation throughout this work. To 
reduce the number of accounted diagrams, we limit our- 
selves to backscattering and forward scattering processes 
only. The Umklapp scattering processes will not be con- 
sidered since they only occur when the distance between 
the parallel FS patches is equal to 7T. We intend to in- 
clude these contributions in a later publication. 

Since the fermionic field theoretical RG is not widely 
known among condensed matter physicists we present 
the method at length and in full detail. We begin by 
presenting our Lagrangian model and the reasons why 
we need to apply the renormalization method in the 
first place for the flat FS. We calculate next the one- 
particle irreducible vertex functions associated with both 
forward and backscattering processes in one-loop order. 
We derive the RG equations for the renormalized cou- 
pling functions and show how they flow to strong cou- 
pling in agreement with the parquet and the other one- 
loop RG approaches. We then move on to calculate the 
self-energy in two-loop order. Using the renormalization 
theory we determine the quasiparticle weight Z along 
the Fermi surface. Subsequently, we add the nonleading 
contributions arising from the nonparquet diagrams. As 
a result of that, we show numerically that the coupling 
functions either flow to large plateau values or approach 
zero. Those nonzero plateau values are however sensi- 
tive to our discretization procedure. Therefore, we can 
neither establish their stability nor completely rule out 
that the plateaus represent some intermediate crossover 
regime associated with the inexact discretization of the 
coupled integral equations. At the end, we conclude 
by discussing how this picture should be altered by the 
renormalization of the Fermi surface whose effects will be 
shown in a subsequent work. 



II. FERMI SURFACE AND LAGRANGIAN 
MODEL 



Our starting point is a strongly interacting 2D electron 
gas in the presence of the flat FS shown schematically in 
Fig.l. In order to keep a closer contact with well-known 
works in one-dimensional physics^, we divide the FS into 
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FIG. 1: The 2D flat Fermi surface with rounded corners. We 
divide it into four regions: two of the solid line type and two 
of the dashed line type. The perpendicular regions do not 
mix in our scheme. 



four regions. We restrict the momenta at the FS to the 
flat parts only. The contributions of the patches perpen- 
dicular to the k x and k y directions do not mix with each 
other in our approach. For convenience, we restrict our- 
selves to the one-electron states labelled by the momenta 
p± = k y and p|| = k x associated with one of the two sets 
of perpendicular patches. The momenta parallel to the 
FS are restricted to the interval —A ^ p|| ^ A, where A 
essentially determines the size of the flat patches. The 
energy dispersion of the single-particle states are given 
by £ a (p) = vf (|rl| — kp) depending only on the mo- 
menta perpendicular to the Fermi surface. The label 
a = ± refers to the flat sectors at p± = ±kp respec- 
tively. Here we take kp — A ^ \pj_\ ^ kp + A, where 
A is some fixed momentum cut-off. The Fermi momen- 
tum kp is not renormalized in our approach and takes its 
nonintcracting value. We also neglect the Fermi velocity 
vf dependence along the FS. This is consistent with the 
fact that the Fermi surface will not be renormalized in 
the present work. 

The Lagrangian L associated with the 2D flat Fermi 
surface is given by 

L= Y ^(a) ff (p)[* 9 *-MbJ-M]Vvk(p) 

p,cr,a— ± 

- y ^ ^2 ISz^asSfi-y - 9iS ay 6ps] 

p,q,k a,/3,<5,7 

x ipl +)g (P + q - k) i>l_ h (k) (q) V( +)ce (p) , 

(2.1) 



where we are following the "g-ology" notation. Here the 
VVi) an d ^(i) are > respectively, the creation and anni- 
hilation operators of particles located at the ± patches. 
The couplings g\ and 52 stand for backscattering and for- 
ward scattering. V is the volume of the two-dimensional 
system which, for convenience, will be set equal to unity 
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FIG. 2: The interacion processes in the model and the cor- 
responding Feynman rules for the vertices. The g\ and gi 
couplings stand for backscattering and forward scattering re- 
spectively. 

from this point on. The graphical representations of the 
corresponding forward and backscattering interactions 
are shown schematically in Fig. 2. In all Feynman dia- 
grams, the single-particle propagators G^K and G^\ will 
be represented by a solid and a dashed line respectively, 
according to their association with the corresponding FS 
patches. 

In setting up a conventional perturbation theory to do 
calculations with this model one immediately encounters, 
in one-loop order, logarithmic divergent particle-particle 
and particle-hole diagramsi^i^ 

n<°> (q , q ± =Q,q\\)= f (p) Gf>j (-p + q) , 

(2.2) 

X (0) (q , q± = -2k F , 9 „) = f G$ (p) G^ (p + q), 

Jp 

(2.3) 

where p = (po,Px,P\\) and f p = f%±^^. 

Using the single-particle propagators from the non- 
interacting part of the Lagrangian we find respectively 



and 




where —A ^ qu ^ A, O = 2vpX is a fixed upper en- 
ergy cutoff and u> is an energy scale which goes to zero 
as we approach the Fermi surface. These infrared di- 
vergent functions will appear infinitely many times if we 




FIG. 3: The Feynman diagrams up to one-loop order for the 
four-point vertex in the backscattering channel. The single 
particle propagators are represented by either solid or dashed 
lines according to their association with the corresponding FS 
patches. 

attempt to do a perturbative calculation of the interact- 
ing one-particle Green's function G(±) (p) or the effective 
two-particle interaction, i.e. the one-particle irreducible 
function T^ 4 ) (pi,p 2 ,P3,P4), to infinite order. This will, 
of course, invalidate any conventional perturbation the- 
ory approach to the problem and, at the same time, will 
make meaningless the direct comparison of our results 
with experiment. While the experiments find in prin- 
ciple finite values for the measured physical quantities, 
our perturbative series expansions are plagued with an 
infinite number of powers of infrared logarithmic singu- 
larities. 

We circumvent this problem following the field theory 
procedure of introducing appropriate counterterms in the 
Lagrangian to render the physical parameters finite in all 
scattering channels 1 ^. In doing so, we eliminate the di- 
vergences order by order in perturbation theory. What 
will become clear soon is that the 2D Fermi surface will 
add important new features in this scheme. The coun- 
terterms are now functions of momenta and vary con- 
tinuously along the FS. In the next section, we begin to 
show how to implement the fermionic field theory RG in 
practice by calculating the renormalized one-particle ir- 
reducible functions associated with both backscattering 
and forward scattering up to one- loop order. 

III. RENORMALIZED COUPLING FUNCTIONS 
AT ONE-LOOP ORDER 

Let us calculate initially the one-particle irreducible 
function T^ 4 ) associated with the backscattering channel 
up to one-loop order. Using appropriate Feynman rules, 
we arrive at 

r{ 4) (px,p 2 ,p 3 ) = -i9i + 2g 1 g 2 U ( -°\p 1 +p 2 ) 
-2g 2 lX (Q) (P2 -p 3 ) + 2g ig2 x W (p2 -pa). (3.1) 

where it will be implicitly assumed from this point on 
that pi = pi +p 2 —p3 by momentum conservation. In ad- 
dition, we follow the same conventions as before concern- 
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FIG. 4: The Feynman diagrams up to one-loop order for the 
four-point vertex in the forward scattering channel. 

ing integrals and energy-momenta representation. The 
diagrams corresponding to these contributions are shown 
schematically in FigEJ 

Let us next perform a similar calculation with the one- 
particle irreducible function I^ 4 ) associated with the for- 
ward scattering channel depicted in Fig0J Using once 
again the Feynman rules associated with our initial Z, 
we obtain 

if } (pi,P2,p 3 ) = -i$a + ff?n<°>(pi + P2 ) 
+.g 2 2 n(°)(p 1 + p 2 ) +g 2 2 x i0) (p 3 -Pi)- (3.2) 

It is clear from those two perturbative series expansions 
that the internal particle-particle and particle-hole lines 
will produce logarithmic singularities if the external mo- 
menta and energies are located at the Fermi surface. To 
deal with this problem and to regularize our perturbation 
expansions, we now invoke the field theory method. 



Being the effective two-particle interaction in the 
backscattering channel, we define the finite one-particle 
irreducible Tf^ such that the renormalized backscatter- 
ing coupling function g 1R (pi\\,P2\\,P3f,u) is given by 

ri 4) (pi,P2,Pa) = -igm (pi,P2,Pa;w) , (3.3) 

where p±± — kp, p 2 ± = —kp, p$± = kp, Pa±_ = —hp and 
the corresponding energy components pio = P20 = 
P30 = 3w/2 and 7540 = —lj/2. Similarly for pu_ = kp 1 
P2± = ~k F , p 3 ± = -k F , p 4± = k F and for p w = p 20 = 
w/2, P30 = 3lu/2 and p^Q = —u>/2, we relate the finite 

one-particle irreducible and the renormalized for- 
ward coupling function g 2 R (j> 1 1| , p 2 m , p 3 1| ; w) to each other 



r 2 (px,P2,P3) = -i92R(Pi,P2,P3;ui) . (3.4) 

We now go back to our initial Lagrangian model and re- 
place the coupling constants g\ and g 2 by the renormal- 
ized coupling functions giR({Pi\\};u) and g 2 R({Pi\\}; u>) 
respectively. If we apply the Feynman rules with our new 
forward and backscattering coupling functions to calcu- 
late both T^ 1 and T^ 4 -* at the Fermi surface we find in- 
stead 



and 



r i 4) = -WlR (Pl||ii>2||,2>3||) + An 2 VF \ I dk \\ [dlR (~ k \\ +Pl\\ +P2\\,k\\,p 3 \\) 92R (pi|| , P2\ , k\\) 



Pi 



+ g2R (k\\,-k\\ +Pl\\ +p 2 ||,Pl|| +P2\\ -P3\\) 91R (pi\\,P2\\,k\\) ] + / dk\\ [2gi R (pi\\,p 2 \\ -p 3 \\ + k \\,k\\) 

Jt> 2 

x 9ir (k\\,p 2 \\,P3\\) - 9ir (Pi\\-.P2\\ ~Pa\\ + k h k \\) 92R (k\\,p 2 \\,p 2 \\ -p 3 \\ + k\\) - g 1R (fc||,p 2 ||,P3||) 
x 92R (pi\\,P2\\ ~Pa\\ + k \\,Pi\\ +p 2 \\ -P3\\)] I In \ , (3.5) 

^ = ~ i 92R (Pl||,P2||,P3||) + 4n 2 VF l J dk \\ [dm +P2\\-,k\\ 1 p 3 \\)92R (Pl\\ , ?2\\ , &||) 

+ 91R ( k \\,~k\\ +p 2 ||,Pl|| +P2|| - P3\\) 91R (Pl\\,P2\\,k\\)] - / dk\\[g2R(pi\\,P3\\ ~ Pl\\ + k \\,P3\\) 

Jt> 3 

x 92R (k\\,P2\\,P3\\ + In , (3.6) 

where all the domains of integration T>i are given explicitly in Appendix A. As it stands, the r^ 4 "* 's are divergent as 
lij — >■ 0, contradicting our earlier assumptions given in Eas. (|3.3|l and 1)3.4(1 . To remedy the situation, we add a new 
term to L 



p,q,k a,/3,(5,7 



[Ag 2 S aS 5 01 - AgiS ai S,3s} ^ {+)s (p + q - k) V,L )7 (k) V(-)/3 (<j) i>(+)a (p) j (3-7) 
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where 

Agin (pi\\,P2\\,Pa\\;v) 



4tt 2 vf 



Pi 



dk \\ [din +Pi\\ +P2\\,k\\,p 3 \\) 92R (pi||,P2||,A||) +9iR (Pi\\,P2\\,k\\ 



x 92R {k\\, -fc|| +p x \\ +p 2 \\,Pi\\ +P2\\ -P3||) ] + / dk \\ [ 2 9iR (Pi||>P2|| -Ps|| + k b k \\) 9m (fe||,P2||,P3||) 



T> 2 



- 91R (Pl\\,P2\\ ~P3\\ + k b k \\) 92R (fc||,P2||,P2|| ~ P3|| + k \\) ~ 92R (Pl||,P2|| ~ P3|| + k \\,Pl\\ + P2\\ ~ P3\\) 
X 9lR (fc||,p 2 ||,P3||)]|ln (j^J , 



(3.8) 



and 



A.927? (Pl\\,P2\\,P3\\;v) = AjT 2 Vp I J dk \\ [9™ {~k\\ +Pl|| +P2\\,k\\,P3\\)g2R (P1||,P2||,*||) +91R (Pl||iP2||j*||) 

x 9m (k\\,-k\\ + +P2||,Pi|| +P2\\ -P3||) ] - / dfc|| [g 2 R (pi||,P3|| ~Px\\ + fe ||>P3||) 

Jt> 3 

X 92R (k\\,P2\\,P3\\ ~Pl\\ + fc ||)]| ln (jjj ■ 



(3.9) 




1R s 



/ -^Ag 2R \^ 



r 



this is the appearance of an energy scale u). All physical 
quantities now depend on this new scale parameter. How- 
ever, the original theory does not know anything about 
u>, i.e., the initial parameters do not depend on the scale. 
This naturally leads us to the renormalization group con- 
ditions for the renormalized coupling functions 



FIG. 5: The additional Feynman rules for the counterterm 
vertices. 



d 

(giR 



m 



(3.10) 



The Agi's are continuous functions of the scale uj and 
of the components of external momenta p\\\ ,P2\\, P3\\ only, 
since their related pn_ components are fixed at the Fermi 
surface. The presence of these new terms in the La- 
grangian results in additional Feynman rules represented 
in FigEJ 

We can adjust the counterterms such that all diver- 
gences are exactly cancelled in our series expansion for 



and T 1 ^' in one-loop order. But the price we pay for 



and 



(g 2 R + A# 2 r) = 0. 

aw 



(3.11) 



The RG flow equations for the renormalized coupling 
functions follow immediately from this. We therefore find 
in one-loop order 



' fa = 4^2^\ J v dk \\[9m\- k \\ +Pi\\ + P2\\,k h p n ) g2 R {Pi h P2ph) 

+ 52fl +Pl\\ +P2||,Pl|| +P2|| -P3||) 9m {Pl\\,P2\\,k\\)] 



d k\\ [2gm (Pi||,P2|| -Ps|| + fc lh fc ||) 9m (k\\,P2\\,P3\\) 
- gm (Pi\\,P2\\ ~P3\\ +k\\,k\\)g 2 R (k\\,p 2 \\,P2\\ - P3\\ + k\\) 

~ g2R {Pl\\;P2\\ ~P3\\ + k \\,Pl\\ +P2\\ ~ P3\\) 9lR (k\\ , P2|| , P3||) ] 



(3.12) 
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and 



d 52fl(Pl||,P2||,P3||) 



dio 



4tt 2 « f { J Vl 



dk\\ [g 2 R (-A)|j +Pi|| +P2||,fc||,P3||) 92R (pi\\,P2\\,k\\) 



9lR (fc||, +Pl\\ + P2\\,Pl\\ +P2|| -P3||) 91R (Pl\\,P2\\,k\\) ] 

dk\\ [g2R (Pl\\,P3\\ ~Pl\\ + k hP3\\) 92R (k\\,P2\\,P3\\ ~ Pl\\ + k \\ 



(3.13) 



. U) 



(2) 




— 5 > 5 + 



counrgnemis 



X' + 



IV. SELF-ENERGY AND QUASIPARTICLE 
WEIGHT UP TO TWO LOOPS 



The field theory RG eliminates the infrared divergences 
order by order in the coupling functions. We do this 
by invoking the usual Feynman rules at each order to 
obtain finite result in terms of the renormalized physical 
quantities. Let us continue with the use of this method to 
obtain the electron self-energy up to second order. Using 
our already modified Feynman rules we have 



FIG. 6: The diagrams for the self-energy calculation up to 
two-loop order. The last diagram of —iT, C ounterterm,s is not 
computed at this stage. It represents the AZ contribution 
which will be added later on. 



These RG equations are exactly equal to the ones ob- 
tained by Zheleznyak et al& using the parquet method. 
In fact, performing one-loop RG calculations is equiva- 
lent to sum the parquet class of diagrams up to infinite 
order. 

Having shown in detail how the RG held theory oper- 
ates in one-loop order we can now proceed to calculate 
with this method the electron self-energy up to two loops. 



jRa 



(P) 



.(i) 

J Ra 



-.(2) 
J Ra 



-* counter terms • 



(4.1) 



We show the diagrams corresponding to each of those 
contributions in Fig[B] The one-loop contributions 

— iTi^l and the first two diagrams of ~i^. countert erms pro- 
duce a shift in kp and they are crucial for the renormal- 
ization of the Fermi velocity and of the Fermi surface 
itselfi 1 ^. Those aspects will be discussed in a subsequent 
paper. We neglect their contributions in the present work 

altogether. Let us therefore concentrate our attention in 

(2) 

the two- loop terms —iTr R(T . We find 



-tS 



(2) / \ 
Ra (P0,P±,P||;w) 



(p - v F (p± - kp)) / dk\\dq\\ [2g 1R (-k + p\\ + q\\ , k\\ , q\\) 



T> 4 



x giR (p\\,q\\,k\\) +2g 2 R(p\\,q\\,-k\\ +p\\ + q\\) g 2 R (k\\,-k\\ + p\\ + q\\,q\\) - giR {p\\, q\\, k\\) 
x g2R {k\\,-k\\ +p\\ +q\\,q\\) - g2R {p\\,q\\,-k\\ +p\\ + q\\) gi R {k\\,-k\\ +p\\ + q\\,p()} 



In 



Q — vp (pj_ — kp) — pa — iS 
vf (p± - kp) -p - iS 



In 



fi — Vf (p_l - k F ) + pa - iS 
vf (p± - kp) +pa - iS 



(4.2) 



r 



where P4 is given in Appendix A. If we take the limits 
u> — > with pj_ = kp or vp (pj_ — k F ) = uj — > at pa = 
our perturbation theory produces an undesirable nonan- 
alyticity at the FS. The way to eliminate this is again to 



add a new term to our Lagrangian of the form 



F 



Pf, 



^\a)a (P) \- W t ~ V F (P-L ~ M] ^(a)a (p) ■ 



p,CT,a— ± 



(4.3) 

With this, the new "noninteracting" Lagrangian now 
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reads 

E ( 1 + F (p\\^)) 

p,cr,a— ± 

X ^(a)c (P) I* 9 * ~ y F (P-L " fcir)] t/>(a)<7 (p) . (4.4) 

Let us define F(pii;a;) = AZ(pii;a;) = — 1. 

Now, in passing, we notice that the main effect of this 
new counterterm is to renormalize the fermion field 



inverse of the renormalized one-particle Green's function 
(p). Let us define T R ^ such that at p± = kp and 
Po = U) 



t rI (Po = u,p± = k F ,p\\) = u. 
Taking into account that 

t rI (p) = z (p\\ ; w ) [Po ~ v f (px ~ k F )] 

v( 2 ) 



k* [po,P±,Pfu) i 



1p(a)a (p) -> Z 2 ; w) ^( a )<r (p) ■ (4.5) 
To determine the function AZ (pu ; coj we consider the it follows from our previous result that 
I 



(4.6) 



(4.7) 



Z(p\\;u) = 1 



1 



32^ JVl 



dk\\dq\\[2gx R (-k\\ + p\\ + q\\, k\\, q\\) gm {p\\,q\\, k\\) + 2g 2 R {p\\,q\\, -k\\ + p\\ + q\\) 



x 92R ( k \\,- k \\ +P\\ +Q\\,Q\\) - 9ir (p\\,Q\\,k\\) g 2R (k\\, -k\\ +p\\ + q\\,q\\) - 92R (p\\, q\\, -k\\ + p\\ +q\\) 

x g 1R (fc||,-fe|| +p||+g||,P||)]ln ■ (4.8) 



This reproduces the result obtained earlier by Kishine 
and Yonemitsui^ within a Wilsonian RG framework. In 
addition to that, making use of Eq. l|4.8|l we obtain the 
following RG equation for the quasiparticle weight 



J 



dZ( Pf 



dui 



= l z {p\\;u) 



where the anomalous dimension 7 is given by 



(4.9) 



l{p\\) 



dk\\dq\\\2g\R {-k\\ +p\\ + q\\ , fc||, q\\) g 1R (py, q\\, k\\) + 2g 2R (p\\, q\\, -k\\ + p\\ + q\\) 



32^v% J Vi 

x g 2R {k\\,-k\\ +p\\ +q\\,q\\) - 9ir (py, q\\, k\\) g 2 R (k\\,-k\\ +p\\ +q\\,q\\) - 92R {p\\ , q\\ , -k\\ +p\\ +q\\) 

x 9iR { k \\,~h +q [l ,p ll )]. (4.10) 

I 



Before proceeding with the full calculation of the RG 
equations for the renormalized coupling functions in two- 
loop order we call attention to the fact that up to now, 
due to the addition of several counterterms, our original 
Lagrangian can be written as 



where g\ R and g 2 b are the bare couplings which are in 
turn defined to be 



L = E Z ^la)a (P) l id t ~ V F (P-JL ~ M] Ao)<7 (p) 
p,cr,a— ± 



i=l 



[g2BS a sSp^ - giB^a-jSps] 



E E 

p.q.k a,/3,<5,7 

x ip\ +)s (p + q - k) ^ ( f _ h (k) i> { _ )0 (q) ip {+)a (p) , 



(4.11) 



JiB 



U z (p*\\) 



(g iR + Ag iR ) 



(4.12) 



Together with the quasiparticle weight, they will render 
the theory finite to all orders in perturbation theory as 
will become clear soon. 
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where the Ag^ r ,lZ° p . with x 



lR{x) 



FIG. 7: The nonparquet diagrams for the renormalized four- 
point vertex in the backscattering channel. 



V. RG EQUATIONS AT TWO-LOOP ORDER 

For convenience, let us rewrite our renormalized La- 
grangian in the form 

L= ^\a)a (P) ~ V F (P-L ~ k F)] ^(a)a (p) 

p,cr,a=± 

- ^2 ISZRSctsSp-y — giRSayfips] 

p.q.k a,/3,<5,7 

x ^l +)s (P + q - k) V(_ )7 (k) tp(-)p (q) ip (+ ) a (p) 

+ AZ ^(a) ( r(P)[^^^(P-L- fc ^)]^(a)<T(P) 
p. a, a— ± 

p.q.k a,/3,<5,7 

x i>\ +)s (P + q k) V>(_ )7 (k) V>(-)/3 (q) ip{+) a (p) • 

(5.1) 

Let us initially consider the diagrams obtained with our 
new Feynman rules generated by this Lagrangian in two- 
loop order. At this point, we only mention that to avoid 
double-counting of diagrams we split up the countert- 
erms into separate blocks. There are several equiva- 
lence relations between diagrams generated by the dif- 
ferent constituent blocks. This will be explained further 
in Appendix B. Taking all this into consideration, the 
contributions due to the counterterms AgiR and Ag 2 R 
cancel exactly all the diagrams which are of the order 
of In 2 (Q/uj) in the vicinity of the Fermi surface. As a 
result, in third-order, the only remaining In (il/w) diver- 
gent contributions are produced by the so-called nonpar- 
quet diagrams. 

We display the nonparquet diagrams for gm in Fig. 7. 
In view of that, to be consistent with the RG condition 
given in Eq. l|3.3[) . we have therefore to redefine Agm us- 
ing the same values for energies and momentum compo- 
nents perpendicular to FS as before. We now find 



-loop 



A<?iH (pin,P2||,P3|i) = A gl R 



A 9 2 lR(°a7 + A Slfl(b) 



2— loop 



a.../, are given explicitly 

in Appendix A. 

We can follow the same procedure to determine Ag 2 R 
in two-loop order. The nonparquet diagrams for g 2 R are 
displayed in Fig. 8. As before, we now get 



2— loop 



A 92R (Pl\\,P2 

+ A 9 2 2r{° c ) P + A 9 2 R {d ) 

. a 2 — loop i a 2 — loo 
+ A 92R(g) + A 9 2 R(h) 

2 — loop 



Agl R l °°P 



A 9 2 2 ruT + A 9 2 2 - l ° OP 



A 9 2 2 r{7 
A 9 2 2 r[T 



a 2-loop 
^ g 2R(f) 
a 2— loop 

A 92R(j) i 



*2R(b) 



(5.3) 



where the Ag 2R ^ p , with x — a...j, are given explicitly 
in Appendix A. 

With the expressions of both Agm and Ag 2 R up to 
two-loop order we can proceed with the calculation of 
the accompanying RG equations for the giRS. Taking 
into account the associated Z factors for the external 
momenta and the RG conditions for the bare coupling 
functions dgiB /duj = together with Eq. (|4.1U|) we obtain 
immediately the two-loop RG equations 



dgiR (P1\\,P2\\,P3\\ 

dm 



2 2~2l(Pj\\)9iR(Pi\\,P2\UP3\\ 

3=1 

dAg iR (pi||,P2||,P3| 



dw 



(5.4) 




G) \ 



FIG. 8: The nonparquet diagrams for the renormalized four- 
point vertex in the forward scattering channel. 
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VI. NUMERICAL RESULTS 

Since it is impossible to solve all these RG equations 
analytically, we have to resort to numerical methods in 
order to estimate how the coupling functions change as 
we vary the scale ui to take the physical system towards 
the FS. Here the basic idea is to discretize the FS contin- 
uum replacing the interval —A pn ^ A by a discrete 
set of 33 points. For convenience, we use uu = S!exp(— I), 
where f2 = 2kp\ with I being our RG step. We will 
choose Q/A = 1. In view of our choice of points for 
the FS, we are only allowed to go up to I « 2.8 in the 
RG flow to avoid the distance w to the FS being smaller 
than the distance between points since the discretization 
procedure no longer applies in this case. 

The choice of the initial conditions at I = in the 
RG equations are, in principle, arbitrary. This is related 
to the fact that one may choose any microscopic model 
to start with in order to study its low-energy proper- 
ties. Since we are most interested in relating our results 
to a repulsive Hubbard model in 2D we initially set the 
couplings equal to a Hubbard on-site repulsive interac- 
tion parameter U in our numerical scheme. As we will 
see shortly, the RG flows are very sensitive to this ini- 
tial value. The physical nature of the resulting state is 
influenced directly by that choice. 

In order to show this, at a first stage, we choose 
9 ib. = 92R = 10 > where g iR = 9iR/nv F . Here, we ap- 
ply the standard fourth-order Runge-Kutta method to 
calculate the flow of the couplings ~g lR and ~g 2 R se ^~ 
consistently. We will present our results for three cases. 
The first one is the one-loop approach. Next, we move on 
to a somewhat intermediate step which includes the con- 
tributions of the self-energy or the Z factor feedback into 
the one-loop RG equations and analyze the correspond- 
ing flows. Then we present the full two- loop order RG 
approach, which includes, in addition to the self-energy 
feedback, the so-called nonparquet diagram contributions 
in the RG equations. We will show that in this final case 
the RG flows are very different from the one-loop case in 
view of the fact that the quasiparticle weight Z becomes 
very strongly suppressed as we approach the FS. 

Later on, we will change the initial condition, choos- 
ing a lower value for the renormalized couplings in or- 
der to contrast it with the previous choice. We take 
~9ir = T)2R = 1- We will see in this case that since 
the quasiparticle weight does not change much from its 
unity initial value, the feedback effect of that factor into 
the RG equations is not very drastic. As a result, the 
two-loop flows will resemble very much the correspond- 
ing one-loop flow for this particular choice of initial in- 
teraction strength. 



A. One-loop RG approach 




FIG. 9: The one-loop RG flow of Sm(pi||,P2||,i'3||) for some 
choices of momenta . 



tially constant (in our case, equal to ten), they acquire a 
distinct momentum dependence as we approach the FS. 
In addition, there are no fixed points in the flow, and 
all coupling values go unequivocally to a strong coupling 
regime. 

This result is consistent with other approaches based 
on a one-loop order perturbative expansion such as the 
so-called parquet method and other RG schemes, which 
predict an instability towards an insulating spin density 
wave state with no sign of nonconvcntional metallic be- 
havior ever to be found in the physical system^. 

However, in all those approaches, the quasiparticle 
weight Z is always equal to unity throughout the cal- 
culation. This should be contrasted with the situation in 
which, due to interactions, the quasiparticle weight might 
approach zero as well. To discuss what happens in this 
case, we move on to the next step, which is the inclusion 




(A.-A-A) 



The one- loop results are depicted in FigslSlandlKJI In 

them, we observe that even though the couplings are ini- 



FIG. 10: The one-loop RG flow of g 2 jj(Pi|| 5 P2|| ,P3||) for some 
choices of momenta . 
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g 1R (0,0,0)-7 


g iR (A,-A-A) 


3. , !• A. .\) / / 


■ 1 i 1 i 1 i 1 i 1 i 



0.0 0.5 1.0 1.5 2.0 2.5 



I 

FIG. 11: The one-loop with the Z factor RG flow of 
9ir(Pi II iP2||iP3||) f° r some choices of momenta. 

of the Z factor feedback into the RG flow equations. 



B. One-loop with the Z factor approach 

For this case, we depict the results in Figsllllandll2l 

In that intermediate approach, which was also discussed 
by Kishine and YonemitsuiS, we observe that the initial 
tendency from one-loop analysis for the couplings to go to 
a strong coupling regime is preserved. However, as the 
suppression of Z becomes prominent the flows seem to 
stop at fixed values. In view of the fact that those values 
are sensitive to our discretization procedure, we can not 
associate those results with the existence of stable fixed 
points. If one increases the number of points necessary 
for discretizing the FS by some factor, the fixed values 



(or the plateaus in the figures) increase approximately 
by the same factor. They should reach a higher bound 
regardless of the number of points in the discretization 
procedure in order to be characterized as true fixed points 
of our RG approach. 

Since this last approach is not yet a complete two-loop 
calculation anyway, we move on to the calculation of the 
full two-loop RG equations for our problem. 



C. The full two-loop approach 

In this final case, we solve Ea. (|5.4() numerically includ- 
ing both the self-energy feedback and the nonparquet di- 
agram contributions for two different choices of initial 
conditions g 1R = g 2R = 10, and g 1R = g 2R = l, 

For the first choice, the results are shown in FigsfEU 
and 1141 In these figures, we observe that the RG flows, 
at least qualitatively, do not change much from the pre- 
vious case. Again, we verify the initial one-loop tendency 
to flow to a strong coupling regime, but as soon as the 
Z effect becomes important the rate of change of the 
couplings reduces abruptly. However, the renormalized 
coupling functions do not saturate as in the previous case 
(the above-mentioned plateaus) . Although the couplings 
seem to reach a plateau regime they now in fact change 
their values continuously at a very slow rate. 

Since these two-loop results are directly associated 
with the effects produced by the quasiparticle weight 
we might as well analyze the flow as we approach the 
FS. In order to do that, we must also calculate Ea. (|4.9() 
self-consistently using the same method described above. 
This equation is solved numerically assuming that for 
I = the quasiparticle weight Z(l = Q) = 1. In doing 
that, we observe numerically a rapid and very anisotropic 
suppression of Z for a moderate coupling regime when we 
approach the FS (see Fig [T3|) . The quasiparticle weight 
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300- 




FIG. 14: The full two-loop RG flow of <?2ii(Pi||>,P2||jP3||) f° r 
some choices of momenta. 



vanishes faster at the center parts of the square FS, i.e. 
in the (jr, 0) or (0, 7r)-directions. This is in qualitative 
agreement with Kishine and Yonemitsu^, although they 
ignore the nonparquet contributions and hence they do 
not perform a full two-loop calculation. 

Besides the rapid suppression of the quasiparticle 
weight Z, the two-loop flow of the coupling functions dis- 
plays other crucial differences in contrast with the one- 
loop case. While in the one-loop order all couplings seem 
to increase indefinitely, this is certainly not the case in 
two-loop. There exists also many couplings which go to 
zero continuously as we approach the FS. In fact, when- 
ever g l (pi||,p 2 ||,P3||) is such that p X || ^ p 2 || + P3||, the 
couplings renormalize to zero and are irrelevant in the 
RG sense. In Fig llfil we show some couplings which be- 




FIG. 16: The vanishing of g 2 ij(Pi|| iVi\\ >i>3||) an d 
5ifl(Pl||)P2||)P3||) f° r some choices of momenta. 



long to this class. 

Now, let us move on to the second choice of initial 
condition for the couplings, i.e. g 1R = g 2 R = 1- Since 
we already know that the distinct two-loop features are 
directly related to what happens to the quasiparticle 
weight, we will limit our RG analysis to the investiga- 
tion of the Z flow in this case. As we see from Fig ll7l 
Z is only mildly different from unity in the RG flow. As 



1.0 



0.9 




0.4 



0.3- 



0.2 



0.0 











v l = 2.0 












1 = 2.16 






I = 2.32 






I = 2.4 




t- 


I = 2.48 




\ 


I = 2.52 




1 = 2.8 
i i i i i i i 



-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

Pi, ia Pi, ,a 




-1.0 -0.5 0.0 0.5 1.0 

P// /A 



FIG. 15: The RG flow for the quasiparticle weight as a func- 
tion of the momentum py along the FS for g 1R (l = 0) = 
9 2 r(1 = 0) = 10. 



FIG. 17: The RG flow for the quasiparticle weight as a func- 
tion of the momentum py along the FS for g 1R (l = 0) = 

g 2R (l = 0) = l. 



12 



a result, the corresponding two- loop flows are very simi- 
lar to the one-loop flow in this case. In order words, as 
expected, for initial weak coupling strengths only, it is a 
reasonable approximation to keep Z = 1 fixed through- 
out the calculation. 

Therefore, we observe that models with different phys- 
ical initial coupling strengths may produce ultimately 
very different physical states. In our numerical estimates, 
there is a clear indication of crossover behavior between 
the regimes defined either by a strongly suppressed quasi- 
particle weight or a Z nearly equal to unity when we vary 
the initial values of the couplings. 



VII. CONCLUSION 

In summary, we showed in full detail a two-loop field 
theoretical RG calculation for a square and fixed two- 
dimensional FS with rounded corners. We neglected 
Umklapp effects at this stage since they only occur when 
the distance between the parallel FS patches is equal to 
7r. We derived explicitly the Feynman rules of the renor- 
malized Lagrangian which regularize all the divergences 
generated by perturbation theory. We did that by con- 
structing appropriate counterterms order by order in per- 
turbation theory. Contrary to what happens in local field 
theory models, those counterterms are now functions of 
momenta varying continuously along the Fermi surface. 
We showed that our method reproduces the results of 
other RG schemes in one-loop order with all the coupling 
functions flowing unequivocally to strong coupling. This 
behavior is typical of the dominant spin-density wave in- 
stability in that regime, which signals antiferromagnetic 
ordering as well as the Mott insulating condition. 

Then, we moved on to the calculation of the self-energy 
correction at two-loop order. Using appropriate assump- 
tions about the form of the inverse of the interacting 
one-particle Green's function at the FS, we were able to 
calculate the quasiparticle weight Z, which was shown to 
vary continuously along the FS. Next, we calculated the 
RG equations for the coupling functions, taking into ac- 
count both the leading order term and the feedback of Z 
into the equations themselves. Although the initial ten- 
dency observed in the one-loop analysis of flowing to a 
strong coupling regime remained, the flows then stopped 
suddenly at fixed plateau values. Despite that, we em- 
phasized that this was not a clear cut evidence of fixed 
points in our approach since those values are sensitive to 
the number of points used in the discretization procedure. 

Finally, we proceeded with the calculation of the full 
two-loop order RG equations for the renormalized cou- 
pling functions. In order to do that, we had to calcu- 
late also the nonlcading divergences which arose from 
the so-called nonparquet diagrams. As we have seen, 
their contribution was not cancelled out by the diagrams 
generated by the counterterms AgiR and Ag 2 R defined 
in one-loop order. That automatically led us to redefine 
the counterterms which from then on acquired also third- 



order contributions. For this final case, we chose two 
different initial conditions to analyze the flows, namely 
9lR = 92R = 10 and g 1R = g 2R = 1. 

For the first choice of initial conditions, the flows were 
qualitatively similar to the ones which included the feed- 
back of Z. They were, as a consequence, very dis- 
similar to the one-loop flow. Afterwards, we numeri- 
cally analyzed the corresponding flow of the quasiparticle 
weight as we approached the FS. We observed that Z was 
anisotropic and strongly suppressed in that limit. We 
pointed out that Z vanishes faster in the (tt, 0) or (0, re- 
directions, i.e the center parts of the square FS. Finally, 
we emphasized that, in view of this strong suppression, 
there were many coupling functions that renormalized to 
zero in contrast to the one-loop case where they appeared 
to increase indefinitely. 

In contrast, for the second choice of initial conditions 
(Sir = 92R = I); % changes very mildly from its orig- 
inal unity value. Consequently, the two-loop flows were 
very similar to the one-loop ones. That is an indication, 
verified numerically, that the physical state depends di- 
rectly on the initial choice of coupling strengths. In view 
of that, we intend to explore further the implications of 
these results. In addition, we must also analyze how the 
renormalization of the Fermi velocity and the FS affect 
our results. This is even more important in view of the re- 
cent results obtained by Afchain et al}- concerning the 
existence of a Luttinger liquid state in two spatial di- 
mensions. As we already noticed from our self-energy at 
one-loop order, the Fermi surface and the Fermi velocity 
are directly affected by interactions. The flows of the FS 
and vf to critical regimes naturally impose important 
new conditions on the flow of the renormalized coupling 
functions. The coupled RG equations for the renormal- 
ized coupling functions, the FS and Fermi velocity should 
all be solved self-consistently. That naturally places ex- 
tra difficulties but one should, nevertheless, deal with this 
problem seriously. We plan to consider those questions 
in a subsequent publication 20 . 
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APPENDIX A 

In this appendix, we will write down the explicit form 
of the nonparquet contributions which are taken into ac- 
count in the Eqs. (|5.2[1 and (|5.3|l . We also give the several 
intervals of integration along the Fermi surface that are 
considered throughout this work which are the following 



V, = 



v* = 



v 4 = 



v, = 



Va = 



-A feii 


<A, 




-A +P2|| 


- fc|| 


-A ^ k\\ 


< A, 




-A sc k\\ 


+ P2|| " 


-P3|| 


-A s$ fc„ 


<A, 




-A sc fc|| 


+ P 3 || - 


"Pill 


-A fc|| 


< A, 




-A q\\ 


< A, 






<A, 




-A ^ -, 


fey +p|| 


+ 9|| 


-A ^ 


< A, 




-A < 9 „ 


<A, 




-A < fe|. 


+ 9|| - 


Pi|| ■ 


-A < fe|. 


+ P211 • 


~P 3 || 



s$ A. 



«S A. 



A. 



< A, 
s$ A. 



A < fe|| < A, 
A < 9|| A, 



-A < k 
-A s? k 



3|| -P2|| < A, 
P 3 || ~P2\\ < A. 



We begin with the nonparquet diagrams shown in 
Fig.Q associated with we get 



32+4 



In 



x / dk\\dq\\g 1R {k\\ + 



n 

to 
Pi I 



•P3||,P2||,fc|| ; 



x giR (pi\\, k \\ +q\\ -P3\\,k\\ +pi\\ -P3\\) 

x 92R (fc||,?||,fe|| +5|| -p 3 ||) , 



(Al) 



a 2 — loop 



1 1 



dk\\dq\\gi R (fc||,<7||,p 3 ||) 

x Slfl ( fc || + Pl|| -P3||,P2||,fc||) 
x ff2fl (Pl||,fc|| +9|| -P3||,9||) , 

(A2) 



a 5i ^7 



1 



1 ^ 

32ttV p m V to 



dk\\dq\\g 2 R (<7||, fc||,P4||) 

F \™ / JT>i 

x 9m (Pi|h fc || +P3|| -Pi||)P3||) 
x 9m (fc|| + <7|| -P4||,P2||,«||) , 

(A3) 



a 2— loop ^ 



In 



x / dk\\dq\igiR (pi||,fc|| +p 3 || -Pi||,p 3 ||) 

X #2i? (fc|| + <?|| -P4||,P2||jfc|| +P 3 || -Pl||) 

X fflil (?|| , fe|| , fe|| + £?|| - P4||) , (A4) 



v 7 = 



-A < fe|, < A, 
-A < 9| , A, 
-A < fe|| +g|| - 

-A SC fe|| +p 3 || 



P4|| < A, 
-Pi|| < A. 



-A < fe|| < A, 

-A < g„ A, 

-A s? fe|| +g|| -p 3 || < A, 

-A s? fe|| + Pn -p 3 || < A. 



In addition to that, we must have of course 



- A < pi], < A, 
-A<p 2 || < A, 
-A < p 3 || < A, 
-A<p 4 |i s$ A. 



Ag 



2 — loop 
lR(e) 



x / dk\\dq\\g 1R (fc|| + pi|j -p3||,P2||, fe||) 
Jt> 8 

X 92R (Pl||)fc|| + 9\\ -P3||,9||) 
x giR (fcj|,g||,fe|| +q\\ -P 3 ||) , 



and finally 



Ag 2-loop _ 



'lR(f) 



dk\\dq\\g 2 R (fey, 9||,P4||) 

>F JT> 7 

X 92R (fe|| +9|| -P4||,P2||,fc|| +P 3 || -Pl|| 

x gm {Pi\\, k \\ +P 3 || -Pi||,P 3 ||) l n 



(A5) 



(A6) 



using the same conventions for the integrals over mo- 
menta components along the Fermi surface as before. 
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Similarly for the nonparquet diagrams shown in Fig. 
associated with F^ we get 



in 

32tt 4 v 2 f \ ... 



x / dk\\dq\\g 1R {p 1 \\,k\\ + q\\ -pi||,fc||) 
Jt> 5 

X QlR ( fc ||)P2||,*|| +P2\\ ~P3\\) 
X 9lR { k \\ +P2\\ ~Ps\\,(l\\,Pi\\) , 



(A7) 



* 2-loop _ 1 



dk\\dq\\g 1R (pi\\, k\\,Pi\\) 



x 9ir +P3|| -P2\\-k\\ +q\\ ~P2\\) 

x 9iR + Q\\ ~ P 2 || ,P2\\ , Q\\) In ^ ~^ > (A8) 



dk\\dq\\g 2 R (fc||,P2||,P3||) 



x S-iij (pi||,fc|| +<Z|| -Pi||,*||) 

x (fc|| + P 2 || -P3||,g||,P4||) ln f^J > ( A9 ) 



a 2~loop 
A 9 2R{d ) = 



16ir 4 Vp \ui 



x / eWciidgH^fl (pi|| , fc|| , fe|| + P3|| - P2||) 

J V 6 

X 9lR {Q\\,k\\ +P3|| -P2|h*|| +9H -P2||) 

x ffiii (k\\ + q\\ -P2||,P2||,9||) , (A10) 



32tt 4 v|, ln (w 



X / *||rfg||52fl (Pl||,fc||,fe|| +P3|| -P2||) 



X fflfi +P3|| -P2||,fc|| + ?|| -P2||) 

x 92R (k + q\\ -p 2 ||,P2||,^||) , (All) 



1 



2 — loop 

2 «(« "32^2 



x / dk\\dq\\g 2 R (piy, fc||, +p 3 || ~P2||) 

x ,<?i/i (A; + q\\ -P2\\,P2\\,q\\) 

X ff2fl (<?|h&|| +P3|| -P2||,P3||) , 



(A12) 



A3; 



2-loop _ ^_ 

2R (s) ~ 32tt 4 i 



dk\\dq\\g 2 R (fc||,p 2 ||,P3||) 



F JT> 5 

P2\\ -P3||,9||,P4|| 



X giR (fc|| 

x 92R (pi\\,k\\ +q\\ -pi||,£?||) hi f^J , (A13) 



a 2 — loop ■ 

* 9 ™W - 32^T j 



dk\\dq\\g 2R (fc||,p 2 ||,P3||) 



X g 2 R +P2|| -P3|h<7|h k \\ + 9\\ ~Pi\\) 



x giR (pi\\,k\\ 



l>lt k il )ln [ " ) . (All) 



Aff 2 2 «;:r = 



1 



167T 4 ? 



x / dfc||dg||ff2fl(pi|j,fc|j,fc|| +p 3 || -P211) 

X ff2fl (fc|| + <?|| -P2||,P2|h*||) 
X g2R (q\\,k\\ +p 3 || -P2||,P3||) , 

and finally 

1 



(A15) 



A 2-loop _ 



J 2R(j) 



16n*v 2 F JVi 



dk\\dq\\g 2 R (fc||,P2||,P3||) 



x g 2R (k\\ +p2\\ -pa\\,q\\,k\\ +q\\ -Pi\\) 

x g 2 R (Pi\\,k\\ +q\\ -Pi\\,q\\) In f^J . (A16) 



APPENDIX B 

In this appendix, we explain briefly how the cancel- 
lation of the In 2 (Q/uj) contributions takes place in two- 
loop order within our scheme. In higher-loop orders, only 
the cancellations of the leading divergence are guaranteed 
by the existing one-loop counterterms. However, we saw 
that in two-loop order there were also new diagrams as- 
sociated with nonlcading logarithmic divergence, which 





FIG. 18: The counterterm diagrams of Agm up to one-loop 
order displayed in terms of constituent blocks. 
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FIG. 19: The counterterm diagrams of Ag2j? up to one-loop 
order displayed in terms of constituent blocks. 





FIG. 20: A two-loop example with Agm mixed with the usual 
interaction vertices. 



Firstly, recall Eq. i|5.1[l , in which the Lagrangian is writ- 
ten in terms of the renormalizcd parameters of the the- 
ory together with the corresponding counterterms. Since 
every term in the interacting part of the Lagrangian gen- 
erates Feynman diagrams in perturbation theory, so do 
the counterterms in our case. 

To illustrate our argument, the counterterms in one- 
loop order for Ag\^ loop and Ag\~j^ oop are now displayed 
graphically in terms of constituent blocks in FigsllSland 
1191 In this way, when we go to higher orders, Feynman di- 
agrams are generated, in which these blocks mix with the 
usual interaction vertices. This will produce the expected 
cancellations of the leading divergence in that order of 
perturbation. As an example, consider the diagrams of 
Figs[2n| and [21 Since the Ag\^ loop and Ag\~ R loop are 
logarithmically divergent on their own, when they are 
inserted in a logarithmically divergent bubble, they pro- 
duce a In 2 (Cl/u) contribution. This cancels exactly the 
leading order divergence produced by conventional per- 
turbation theory. However, in order to avoid double- 
counting of these counterterms contributions, we need to 
take into account the various symmetry relations which 
exist among diagrams constructed with some of the coun- 
terterms blocks. We illustrate one of these symmetry 
relations in Fig[23 Following this, only the single log- 
arithmic divergences from the nonparquet diagrams are 
left out of these cancellation processes and they are used 
to redefine Agm and Ag2R in two-loop order. 



FIG. 21: A two-loop example with Ag2R mixed with the usual 
interaction vertices. 



were not accounted by the one-loop order counterterms. 
This led us to redefine Agm and Ag2R in two-loop order. 
This is consistent with the fact that we do not know, a 
priori, the exact expression for the counterterms and, for 
this reason, they must be calculated order by order in 
perturbation theory. 



FIG. 22: A symmetry relation example of diagrams generated 
by different constituent blocks. 
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